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ABSTRACT 



The transmission of light through a planetary atmosphere can be studied as a func- 
tion of altitude and wavelength using stellar or solar occultations, giving often 
unique constraints on the atmospheric composition. For exoplanets, a transit yields 
a limb-integrated, wavelength-dependent transmission spectrum of an atmosphere. 
When scattering haze and/or cloud particles are present in the planetary atmo- 
sphere, the amount of transmitted flux not only depends on the total optical thick- 
ness of the slant light path that is probed, but also on the amount of forward- 
scattering by the scattering particles. Here, we present results of calculations with 
a three-dimensional Monte Carlo code that simulates the transmitted flux during 
occultations or transits. For isotropically scattering particles, like gas molecules, 
the transmitted flux appears to be well-described by the total atmospheric opti- 
cal thickness. Strongly forward- scattering particles, however, such as commonly 
found in atmospheres of Solar System planets, can increase the transmitted flux 
significantly. For exoplanets, such added flux can decrease the apparent radius of 
the planet by several scale heights, which is comparable to predicted and measured 
features in exoplanet transit spectra. We performed detailed calculations for Titan's 
atmosphere between 2.0 and 2.8 /zm and show that haze and gas abundances will 
be underestimated by about 8% if forward-scattering is ignored in the retrievals. 
At shorter wavelengths, errors in the gas and haze abundances and in the spectral 
slope of the haze particles can be several tens of percent, also for other Solar Sys- 
tem planetary atmospheres. We also find that the contribution of forward-scattering 
can be fairly well described by modelling the atmosphere as a plane-parallel slab. 
This potentially reduces the need for a full three-dimensional Monte Carlo code 
for calculating transmission spectra of atmospheres that contain forward-scattering 
particles. 



Keywords: atmospheres, composition ; extra-solar planets ; Titan, atmosphere ; ra- 
diative transfer 



3 



1 Introduction 



Planetary atmospheres can be studied remotely using spectroscopy of the light that 
is reflected or the radiation that is emitted by the planet. Another method for study- 
ing planetary atmospheres is to measure the light of a bright source, such as a star 
or the Sun, as it is attenuated by (part of) the atmosphere. For Solar System objects, 
such transmission measurements are performed during stellar or solar occultations, 
when the planetary limb is in between the light source and the spacecraft or tele- 
scope. One advantage of transmission measurements is that the light source is often 
very bright, enabling high signal-to-noise measurements. Another advantage is that 
transmission measurements allow sampling of long path lengths through the at- 
mosphere, since the light travels through the curved limb of the planet. Long path 
lengths increase the sensitivity to e.g. trace gases with very small concentrations. 
Furthermore, altitudes in an atmosphere can be probed with a high vertical resolu- 
tion, as only a very limited altitude range is probed by a single observation com- 
pared to e.g. on-disc observations, and many altitudes can be sampled due to the 
high signal-to-noise. Also, occultation measurements rely on relative changes in the 
measured signal, and not on absolute signal levels, and are hence 'self-calibrating' 
and less dependent on instrument drifts. As a result of these advantages, transmis- 
sion measurements are very suitable for retrieving vertical profiles of trace gases. 
Disadvantages of solar and/or stellar occultations are that they require very spe- 
cial observing geometries, making it not trivial to target specific places on the 
planet. Furthermore, only a single location on the planet can be studied simulta- 
neously, and the long path lengths prevent probing the lower, thicker layers of a 
planetary atmosphere. For solar occultations, the atmosphere can naturally only 
be studied at local twilight conditions. In recent years there have been many oc- 
cultation measurements of Solar System planets. Besides the Ea rth, occultation 



2008; 


Fedorova et al., 


2008: 


Krasnopol'Skii et al.. 


1989; 



Belvaev et al. 



Forget et al. 



2012), Mars (e.g. Blamont et al, 1989 



20091). Jupiter (IFormisano et all 120031) . 



Saturn dKim et all l2012h. Titan dBellucci et all lipOft iKim et all 1201 1|). and Pluto 



(Hub bard et all 119881 : lElliot et alll2003|) . for which more are planned (|Stern et al. . 
2008). Besides measuring the atmospheric transmission, occultations can also be 
used to measure the atmospheric refraction, which we will not consider here. 



Transmission measurements are one of the most valuable probes of the atmospheres 
of extrasolar planets (exoplanets). The transmission of an exoplanetary atmosphere 
can be measured when the orbit of the exoplanet with respect to the observer is 
such that the planet crosses the disc of its host star. By accurately measuring the 
wavelength-dependence of the decrease of star light during s uch a transit, a trans- 
mission spect r um of the planetary atmosph ere can be derived (ISeager and SasselovL 
2000l : iBrownl l200ll : iHubbard et all l200ll) . although only the day-night boundary 
can be probed. These transit measurements have revolutionised our knowledge of 
exoplanet atmospheres and are now widely used to derive atmospheric properties of 
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a number of transiting exoplanets. There are now several planets for which molec- 
ular absorption fea t ures have been identified through measuring their transits (e.g. 
Charbonneau et al. . 2002 ; Tinetti et al. . 2007 : Snellen et al. . 2010 ; Seager and DemingL 
2010) , as well as molecular and haze Ray leigh scattering (e.g. lLecavelier Pes Etangs et all 



2008 



Pont et all 120081 : ISing et all 1201 II). Also, there is a n indication of high wind 



speeds on the planet HD 209458b (ISnellen et all |2010|) . The smallest planet for 
which there have been significa nt transit measurements is the 'super-Earth' GJ 



1214b (Char bonneau et all 12009). T ransit measurements allow the d etermination 



of the bulk composition of this planet (Mil ler-Ricci and FortnevLl2010() . but at present 
it is not clear whether it i s has a water vapour atmosphere, or a hydrogen atmos phere 
with clouds dBean et al.l . [20 1 d.l20 1 ll : Icroll et all EoTTl : Ide Mooij et alll2012h . Un- 
fortunately, atmospheres with high molecular weight, like Earth, are hard to charac- 
terise in this way, as the spectral features are much reduced compared to hydrogen 
atmospheres. This is because the scale height of the atmosphere scales with the 
molecular weight. 

The transmission of a planetary atmosphere is mainly determined by the extinc- 
tion or total optical thickness of gases and particles, like haze and cloud particles, 
measured along the line of sight. The total optical thickness r t can be split into an 
absorbing part, r a , and a scattering part, r s , using the single-scattering albedo u of 
the mixture of gases and particles in the layer: 



r t = ur t + (1 - u)r t = r s + r a 



(1) 



The optical thicknesses r t , r a , and r s and the single scattering albedo u are usually 
wavelength dependent. 

If the transmission T, i.e. the ratio of transmitted to incident flux, is purely de- 
termined by the extinction along the line of sight, it is simply given by the Beer- 
Lambert law: 

T = e~ Tt (2) 



If the scattering optical thickness r s along the line of sight is larger than zero, how- 
ever, a fraction of the light that has been scattered out of the incident beam of light, 
will be scattered in exactly the forward direction after a single scattering or after 
multiple scattering events. This means that more light reaches the observer than ex- 
pected from the extinction optical thickness of the atmosphere alone. In fact, large 
particles scatter light predominantly in the forward direction because of diffraction 



(see e.g. Ha nsen and Travisl |1974|) . possibly making the forward-scattered light 
contribute significantly to the total transmitted flux. Indeed, non- absorbing, per- 
fectly forward-scattering particles would appear completely transparent in trans- 
mission, despite a non-zero optical thickness. 

How a particle scatters the incident flux as a function of the scattering angle 6 
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is described by the phase function P(6). Figure 1 shows the single- scattering 
albedo Co and the phase function in the forward direction, P(0), of some parti- 
cles found in atmospheres of Solar System planets. The phase functions have been 
normalised such that their integral over 4tt steradians is unity. We show measured 
as well as calculated values of P(0). The properties of Titan aerosol and mar- 
tian dust particles ar e derived from measurements by iTomasko et al.l (|2008) and 
Tomasko et al.l (119991). respectively. The lines in Fig. 1 are computed using Mie- 
theory (|de Rooij and van der Stap . 1984 ). Note that the strength of the forward- 
scattering part of the phase function is due to diffraction and as such depends 



most ly on the size of the particles, and less on their shape (e.g. iMishchenko et al. 



1996). Hence, Mie scattering should be accurate for the purpose of estimating the 
forward- scattering peak. For the Ven us cloud particles, we used the optical con- 
stants of Palmer and Williamsl (11975b and the size distributions for the two most 
dominant size modes from iGrinspoon et al.l (I1993|) . For the calculations for mar- 
tian dust and ice clo ud particles we used th e optical parameters and size distribu- 
tions as described by IKleinbohl et al.1 (|2009l) . It can be seen that especially at short 
wavelengths P(0), and hence the forward-scattering contribution to the transmis- 
sion, can potentially be large, even though all these different particles are relatively 
small (with effective radii smaller than 2 fim). The forward-scattering part of the 
phase function of larger particles, such as liquid water cloud particles found on 
Earth, can be much larger (see e.g. Fig. 11.13.6 in Mishch enko et al.U2006|) . 
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Fig. 1. Single-scattering albedos uj and phase functions in the forward direction, -P(O), of 
particles present in the atmospheres of Venus, Mars and Titan (see the text for details). 

Analysing transmission measurements using Eq. 2 and thus ignoring the (poten- 
tial) contribution of forward- scattered flux to the transmission, which is how these 
measurements are usually analysed, will result in an error in the retrieved optical 
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thickness of the planetary atmosphere. In this paper, we explore the contribution 
of the forward- scattered flux to the transmission signal using a three-dimensional 
Monte Carlo model, which simulates scattered light in a spherical, stratified atmo- 
sphere. Besides providing general results, we will also present more detailed results 
for Titan's atmosphere, to serve as a practical example. Finally, we draw conclu- 
sions regarding possible errors in retrievals from transmission measurements when 
scattering particles are present. 



2 Numerical model 



Our model is a Monte Carlo model, whose core is based on subroutines of the 



mc-u nik code by Andreas Macke (pers. comm.: iMacke et all 1 19991 : ICahalan et al. 



2005). In a Monte Carlo model, the light that is incident on the planetary atmo- 
sphere is described by a large number of separate photons, the paths of which are 
followed through the atmosphere. In our model, photons are fired through the at- 
mosphere from one side and in one direction, and their locations and directions 
when they leave the atmosphere (provided they have not been absorbed) are stored. 
We ignore refraction and polarisation. Refraction will generally increase the path 
length through t he atmosphere, b ut it s effect at low pressu res has been found to 



be small by Hubbard et al. ( 2001 ) and Bellucci et al. ( 2009|) . A model atmosphere 



is assumed to be spherically symmetric and is described by an arbitrary number 
of shells. In our calculations, the thickness of each atmospheric shell is roughly a 
sixth of the atmospheric scale height. We did not try to find the optimum for this 
thickness, but halving the number of layers did not significantly change our scat- 
tering results. Also considering our simulations with a single layer, which show 
that the spherical geometry is not important to first order for the scattering con- 
tribution (see Sections 3 and 4), we do not expect large errors in the scattering 
contribution due to our finite number of layers (of the order of 50). When analysing 
real measurements, more layers might be needed, depening on the required model 
accuracy. Across each atmospheric shell, the temperature, pressure, gas and parti- 
cle concentrations and their optical properties are assumed to be constant. In the 
results presented below, we always have only one type of scatterer present in the 
atmosphere. The atmosphere is bounded below by a black surface at a pressure of 
5 bar. Generally, the slant optical thickness of a light path crossing the 5 bar level 
is of the order of 100 or more and hence the light emanating from this level has a 
negligible contribution to the transmitted signal. 

When a photon enters an atmospheric shell (either the outer or an inner shell, and 
either from the outside or the inside), the free path length /* of the photon in the 
shell is calculated using: 
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where r* follows from 



(4) 



with £ a random number between and 1 . Within each shell in our atmosphere, 
dr t /dl is constant. 

If the free path length and the direction of the photon are such that the photon passes 
through the shell without being intercepted by a gas molecule or a particle, we stop 
the photon where it would leave the inner or outer boundary of the atmospheric 
shell. If the photon is leaving the outer atmospheric shell towards space, its location 
and direction of propagation are stored and a new photon is fired. If the photon is 
about to enter another atmospheric shell, we calculate the photon's free path length 
in the new shell and let it continue its path in the new shell. 

If the free path length and the direction of a photon in an atmospheric shell are such 
that the photon is intercepted by a gas molecule or a particle in the same shell, we 
calculate whether the photon is absorbed or scattered, based on the single-scattering 
albedo u> of the mixture of gas molecules and particles in the shell. If the photon is 
absorbed, it is lost and a new photon is fired. If the photon is scattered, its direction 
changes according to the phase function P of the mixture of gas molecules and 
particles in the shell. The photon's scattering angle 6*, and hence its new direction, 
is determined by the cumulative distribution of the phase function as follows: 

e* 7T 

J P(G*)de* = i J P(Q)dQ, (5) 
o o 



with £ a random number between and 1 . For this new direction, a new free path 
length I* is calculated and the photon is sent on its way. 

A photon reaching the surface below the atmosphere is assumed to be absorbed. 
So, all photons are either absorbed or leave the top of the atmosphere. Figure 2 
shows a sketch of the geometries involved for photons fired at a tangent altitude 
h. The figure includes sample trajectories of photons calculated using our Monte 
Carlo model. For exoplanets, illumination from all values of h, reaching from the 
bottom to the top of the atmosphere, has to be taken into account. In our model, 
we find that the contribution of photons exiting at a different tangent altitude from 
where the photons are inserted, in exactly the right direction, is negligible. 

We have tested our Monte Carlo code using an atmosphere consisting of a plane- 
parallel slab, for which several accurate codes for calculating scattered fluxes exist. 
We compared our calculations of the total transm itted flux and the flux t ransm itted 



in a specific direction for isotropic scattering with Sh ettle and Weinman! (11970 ) and 



de Haan et al.l (| 19871) (ignoring polarisation), and for non-isotropic scattering with 



8 



Fig. 2. Illustration of the geometries. The illumination source is on the left side of the 
planet, and the observer on the right side. Photons are fired in a single direction on the 
left, indicated by the arrow on the left. The observer measures the photons that leave the 
atmosphere in the direction of the arrow on the right. The solid lines indicate the top and 
bottom of the atmosphere in the plane of the emission and observer and the tangent altitude 
of the observation is denoted by h. Apparent penetration of the rays at the bottom of the 
atmosphere is caused by the rays moving in the dimension perpendicular to the figure. 



de Haan et al.1 (|l987r) and iHovenier and de Haanl (|1985f ). In all cases, the results 



agree to within a percent. We have also tested our code assuming a spherical shell 
model atmosphere containing only purely absorbing particles. For such an atmo- 
sphere, the number of photons transmitted through the atmosphere, divided by the 
total number of fired photons, should simply equal the transmission T as described 
by Eq. 2, which indeed it does to well within a percent. Hence, our sampling of the 
mean free paths is also accurate in the three-dimensional model. 

For every tangent altitude h, we fire a number of photons that is equal to 10 5 divided 
by the transmission derived from the total optical thickness along the line-of-sight 
(Eq. 2), to ensure that enough photons survive their trip through the atmosphere to 
have relatively accurate results. Afterwards, we normalize the number of surviving 
photons to the number of incident photons for each tangent altitude h. We impose 
a maximum number of photons of 10 8 for large optical thicknesses, to keep the 
computation time manageable. This means that, for large optical depths, less than 
10 5 photons are transmitted and the accuracy of the result decreases. 



3 Scattering results and application to exoplanets 



We have performed calculations of the flux that is transmitted through a plane- 
tary atmosphere and that is detected at the other side, with the transmitted photons 
travelling in the same direction as the incident ones. In this section, we consider 
three model atmospheres that differ in the type of scattering particles: Rayleigh 
scattering particles, and two types of forward- scattering particles, with different 
forward- scattering strengths. As briefly discussed in the introduction, one would 
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qualitatively expect more photons to reach the observer if scattering is included, 
with more forward-scattering particles giving a larger added signal. In each model 
atmosphere, the particles are uniformly mixed throughout the atmosphere in hy- 
drostatic equilibrium and they are the only opacity source. The three model at- 
mospheres have the same extinction optical thickness, and are based on the 'hot 
Jupiter' HD 18973 3b, for which Rayleigh-scattering particles have been identified 



in the atmosphere (|Pont et all 120081 : ISing et all |2011|) . Although we present here 



only the results with parameters from one specific atmosphere, it will be clear later 
that these results are also valid for other planets. A uniform temperature of 1000 K 
and a mean molecular weight of 2.2 gram/mol is assumed, with a gravitational 
acceleration of 18 m/s 2 . 



The f orward-scattering particles have Henyey-Greenstein phase functions (IHenyey and Greenstein , 
194l|) . given by 



P(6) 



1 



4tt (1 + g 2 -2c/cos6) 3 / 2 ' 



(6) 



with g the asymmetry parameter. For one type of particles, we chose g = 0.9, which 
represents a case of moderate forward scattering, with P(0) = 15.1 (cf. Fig. 1). The 
other type of particles is more forward-scattering, with g = 0.98, corresponding to 
P(0) = 393.8. All particles are assumed to be non- absorbing, with single- scattering 
albedos equal to unity. 

Figure 3 shows the transmitted flux for different tangent heights h (see Fig. 2), plot- 
ted as a function of the total optical thickness along the line of sight for the three 
model atmospheres. In the case of Rayleigh scattering particles, only very little 
extra flux is transmitted as a result of the forward- scattering and the total transmit- 
ted flux can be well described by Eq. 2. Indeed, adding the forward- scattered flux 
to the transmission changes the transmitted flux by less than a percent (relatively 
speaking, so not in percentage points; i.e. a 1% change of a transmission of 20% 
gives 19.8% here, not 19%) for optical thicknesses less than unity, and by a few 
percent for optical thicknesses around 10. In the cases with the forward- scattering 
particles, however, the forward- scattering adds significantly to the transmitted flux 
as expected based on extinction alone (see Eq. 2 and the line for g = 0.0 in Fig. 3). 
Especially at large optical thicknesses, these scattering particles can yield a trans- 
mitted flux that is several times the flux that would be transmitted if the particles 
were purely absorbing. For g = 0.98 this is most apparent, but also for g = 0.9 this 
is the case for very low transmission values, which are not well visible in the fig- 
ure. These results confirm our qualitative expectations, as noted in the introduction, 
but now the magnitude of the effect is made apparent. In terms of absolute number 
of photons, optical thicknesses around unity give the largest differences between 
forward-scattering and purely absorbing particles. From our code, we also know 
the number of times a photon is scattered before it exits the atmosphere at the top 
or the bottom. As expected, for small optical thicknesses, the additional scattering 
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contribution comes mostly from photons that are scattered only once, whereas for 
large optical thicknesses multiple scattering is more important. 




Fig. 3. The transmitted flux as a function of the slant optical thickness through the atmo- 
sphere for isotropically scattering particles (g = 0.0) and two types of Henyey-Greenstein 
scattering particles (g = 0.9 and g = 0.98). The line for isotropically scattering particles 
is practically indistinguishable from a line (not shown) representing transmission without 
added forward-scattering scattering (Eq. 2). 

To simulate the transmission for transiting exoplanets, the transmitted flux must be 
integrated along the planetary limb that is in front of the star. By measuring the 
transmitted flux at a given waveleng th, one can determine the apparent radius of 
the planet at that wavelength (e.g. iBrownl EoOlh . Our calculations show that the 



apparent radius of exoplanet HD 189733b with an atmosphere filled with forward- 
scattering particles, decreases 0.7% and 3.8% for g = 0.9 and g = 0.98, respec- 
tively. This corresponds to roughly 1 and 4 times the atmospheric scale height and a 
decrease in transit depth of 0.03 and 0.15 percentage points (from an initial 2% tran- 
sit depth). Such differences are easily comparable to the depth of spect ral features 



in the transit spectrum, and to differences in measured transit depth (e.g. lSing et al. 



201 lb . For other exoplanets, the number in terms of scale height will be similar, as 



we found by changing gravity. 

We also compared ou r Monte Carlo results with calculations of light passing through 



a plane-parallel slab (Ide Haan et all |1987|) with the same scattering particles and 



the same total optical thickness as the light path in the three-dimensional, spheri- 
cal atmosphere. For all tangent heights and corresponding optical thicknesses, the 
relative difference between the flux transmitted through the three-dimensional at- 
mosphere and that transmitted through the plane-parallel slab calculations is less 
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than 10%, meaning that plane-parallel slab calculations (e.g. Hu bbard et all 1200 1|) 
can describe the effects of forward- scattering particles fairly well. This is perhaps 
not surprising, since by far the largest optical thickness along the line of sight is 
located at a relatively small region in the atmosphere, around the tangent point, for 
each (almost) parallel beam that is passing through the atmosphere, from the star 
to the observer. 

The outcome of this last test suggests that the exact geometry is not important to 
first order, which means that Fig. 3 is valid for all different kinds of atmospheres, if 
there is no absorption. Different atmospheres will then only differ in the amount of 
slant optical thickness as a function of height and wavelength. The case when there 
is also absorption by e.g. gases present is illustrated in the next section. 



4 Application to Titan 



To further explore the effect that scattering haze and/or cloud particles have on 
transmitted fluxes and the atmospheric information derived from transmission mea- 
surements, we take Titan's atmosphere as example, for which a sequence of near- 



infra red solar occultation measurements by Cassini VIMS were analysed by lBellucci et al 
(|2009|) . The presence of gas absorption bands makes this case more complex. Again, 
we expect scattering from the haze to increase the amount of observed photons. In 
gas absorption bands this effect is expected to be less important, but it is not directly 
apparent what the magnitude of the effect is and w hat effect the scattering will have 
on retrieved gas abundances. iBellucci et all (|2009 ) derive the optical thickness and 
spectral slope of Titan's haze and gas abundances, all as functions of altitude, from 
these observations at a single location. Here, we do not attempt to fit the VIMS data, 
but instead simulate transmission spectra using our Monte Carlo code and fully in- 
cluding forward- scattering particles. From these simulated spe c tra, w e then derive 
gas and haze properties ignoring scattering, like lBellucci et al.l (|2009l) do. We then 
compare the retrieved atmospheric parameters with the 'real' parameters that we 
put into our code to quantify potential errors in retrievals when forward- scattering 
is ignored. 

For Titan's atmosphere we calculated tr ansmission spectra b etween 2.0 and 2.8 fim 



with a spect ral resolution of 16.6 nm (IBrown et al.l . |2004|) using the correlated 



k method (ILacis and OinasL |1991|) and CH 4 line parameters from HITRAN 2008 
(|Rothman et all |2009|) . Opacity distribution tables are calculated for the ranges of 
temperatures and pressures relevant for Titan's atmosphere. Since we are investi- 
gating the effects of ignoring forward- scattering and because we are not fitting real 
measurements, we don't have to be concerned about potential errors in the gas ab- 
sorption coefficients due to the low temperatures in Titan's atmosphere; we use the 
same gas absorption properties in the calculations that include forward- scattering as 
in the calculations that ignore this scattering. Similarly, we are not concerned about 
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small errors expected from the correlated- A; approximation, especially since by far 
the most opacity along the light path crossing the atmosphere is located in a sin- 
gle atmospheric layer, whose altitude is the tangent height. Potentially uncorrelated 
/c-distributions at different altitudes will therefore have little effect. We assumed a 



pressure-temperature profile from Vinatie r et al] (120061) and a uniform CH 4 Volume 
Mixing Ratio (VMR) of 1 .4% ([Niemann et all 120051). The haze op tical thickness as 



a function of altitude at 2 /im was taken from lBellucci et a l. (2009 |) and the s pectra l 



spt 

slope of the haze is taken to be that of the upper haze of iTomasko et al.l (|2008|) . 
The single-scattering albe do of the haze particle s as a function of wavelength was 
linearly extrapolated from lTomasko et al.l (I2008|) and the pha se functions wer e lin- 
early interpolated from the phase functions tabulated by lTomasko et al.1 (|2008) (see 
Fig. 1). For our Mon te Carlo calculation s we took into account the distance of the 
spacecraft to Titan (Bellucci et alll2009|) as well as the VIMS telescope dimensions 
(|Brown et all 12004) to collect all of the 'observed' photons. How the photons are 
imaged after hitting the primary mirror is not taken into account, but this is ex- 
pected to have very little influence, because by far the largest contribution of the 
signal comes from the direction of the Sun. 

Figure 4 shows calculated transmission spectra with and without including the 
forward- scattering contribution for different obser ver's tangent a l titude s. This part 
of the spectra shows a CH 4 absorption feature and lBellucci et al.l (I2009f) derive the 
CH 4 VMR from this part of the spectrum (cf. their Fig. 7). From our figure, it is 
clear that the contribution of forward-scattering by the haze particles is largest in the 
continuum, where there is little gas absorption, as expected. The added contribution 
of forward- scattering to the continuum flux makes the absorption features deeper 
with respect to the continuum. It also slightly flattens the slope of the continuum 
flux, because the haze particles are more forward- scattering at smaller wavelengths 
than at larger wavelengths. 

We derived the atmospheric parameters by fitting the transmission spectra that we 
calculated using the Monte Carlo code, which includes scattering. These spectra 
(dashed lines in Fig. 4) are thus thought of as the 'real' measurements that need fit- 
ting. For the fi t ting, w e ignored the scattering, and used just Eq. 2, similar to what 
Bell ucci et al. (2009) do. The retrieved atmospheric parameters can then be com- 
pared to the 'real' state of the atmosphere, as input in our scattering model. We first 
scaled the amount of haze particles and the spectral slope (the Angstrom exponent) 
of the haze to exactly match the beginning (2 /mi) and end (2.8 /im) of the trans- 
mission spectra. Compared to the 'real' values that we used in the calculation of 
the transmission spectra, ignoring the forward-scattering in the fit scales down the 
density of haze particles by typically 8% and the Angstrom exponent by 3%. Then, 
we used a least-squares fit to determine the CH 4 VMR. Figure 6 shows the retrieved 
CH 4 VMR as a function of altitude together with the 'real' VMR. The noise on the 
retrieved values is the result of the limited number of photons used in the Monte 
Carlo calculations. The decrease of the error in the retrieved VMR with altitude 
above about 200 km, corresponds with the decrease of the difference between the 
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Fig. 4. Calculated transmission spectra with (dashed lines) and without (solid lines) for- 
ward-scattering for Titan's atmosphere, for four different observed tangent altitudes (see 
Fig. 2.) 

transmission spectra in Fig. 4. Figure 6 shows that in this case, ignoring forward- 
scattering when retrieving the gas abundance can lead to an underestimation of the 
gas VMR of typically 8%. 



The reason for the underestimation is that ignoring the added continuum flux that 
is due to forward-scattering leads to a lower retrieved haze particle density (since 
more flux is transmitted in the continuum). It is not straightforward to argue what 
effect this will have on the depth of the absorption band compared to the 'real' spec- 
trum, since the scattering contribution to the transmitted signal is not related to the 
optical depth in a simple way as the non-scattering signal (Eq. 2) is. The outcome 
of lowering the haze opacity in the non-scattering case is illustrated in Fig. 5. The 
dotted line shows the fit when the haze is scaled down, but the methane abundance 
is kept identical. The apparent increase in depth of the absorption band between the 
solid (no forward scattering) and dashed (with forward scattering) spectra corre- 
sponds simply to an almost uniform lowering of the optical depth with wavelength, 
giving rise to a non-uniform increase in transmission according to Eq. 2. However, 
whether this increase of band depth is more or less than the increase of band depth 
obtained from the scattering contribution in the dashed spectrum is not immediately 
clear. It can be seen from Fig. 5 that in our calculations, the absorption band is too 
deep for the non-scattering case (dotted line) and the CH 4 abundance should be 
decreased to fit the dashed line. Surprisingly perhaps, ignoring forward-scattering 
thus leads to an un Jerestimation of both the gas and the haz e abundances if mea- 
surements are fitted. It is therefore all the more puzzling that Bellucci et all (|2009r) 
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find a CH 4 abundance that is significantly lar ger than that measure d by the Huygens 
probe ( Niemann et al. . 2005 ) below 200 km. iBellucci et al. ( 20091) suggest various 
possible explanations for this effect, such as extra absorption within the haze itself. 



0.50 



119 km 




2.6 2.7 
Wavelength (micron) 



2.8 



Fig. 5. A zoom of a spectrum with dashed and solid lines as in Fig. 4. The dotted line now 
shows the absorption-only fit to the spectrum when the haze opacity and its spectral slope 
is decreased, but the methane abundance is kept identical. 

As we did for the exoplanet calculation (see Section 3), we compared the trans- 
mission spectra as calculated with our Monte Carlo code and assuming a spherical 
model atmosphere with spectra calculated assuming a plane-parallel slab model at- 
mosphere with the same particle scattering properties, weighted by the slant optical 
thickness of each atmospheric layer, and the same optical thickness (as measured 
along a slant light path for a given tangent altitude). The differences between the 
spherical atmosphere spectra and the plane-parallel slab spectra are small: of the 
order of a few percent in the transmitted flux. 



5 Conclusions and discussion 



We have developed a three-dimensional, spherical Monte Carlo code to simulate 
the added contribution of scattered light to the flux that is transmitted through a 
planetary atmosphere during stellar or solar occultations and exoplanet transits. 
Our calculations with parameters of the exoplanet HD 189733b show that Rayleigh 
scattering does not add significantly to the transmitted flux that one would expect 
when the transmission depended solely on the extinction (absorption plus scatter- 
ing) optical thickness of the atmosphere. However, forward-scattering particles can 
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Fig. 6. The CH4 VMR as retrieved while ignoring forward-scattering by the haze particles 
(solid line). The vertical dashed line indicates the 'real' VMR value that we used as input 
for our spectral calculations, which included forward-scattering. 

contribute significantly to the transmitted flux. For (slant) optical thicknesses larger 
than unity, the added contribution of forward- scattered flux can be as large or larger 
than the flux expected from the extinction alone. For exoplanet transits, this means 
that the magnitude of a transit of a planet with forward- scattering particles in its 
atmosphere can be several percent smaller than predicted when these particles are 
assumed to be fully absorbing. Such a difference in transit depth is easily com- 
parable to the magnitude of measured and predicted spectral features in a transit 
spectrum. 

With our Monte Carlo code, we performed detailed calculations of occultation 
spectra of Titan's atmosphere between 2.0 and 2.8 /im and show that the haze 
abundance and trace gas mixing ratios are underestimated by ~ 8% if the forward- 
scattering contribution to the transmitted flux is ignored in the retrieval algorithm. 
Also in such a case, the derived spectral slope of the haze particles is underes- 
timated by a few percent. Figure 1 shows that Titan's haze particles are fairly 
forward- scattering between 2 and 3fim, but that cloud and haze particles can be 
at least an order of magnitude more forward-scattering at shorter wavelengths. This 
means that at shorter wavelengths, cloud, haze or gas abundances could be under- 
estimated by several tens of percent if forward- scattering is ignored. Figure 1 also 
shows that such strong forward- scattering is not uncommon in other atmospheres 
in the Solar System. Also the change of forward-scattering with wavelength can be 
much stronger at shorter wavelength, resulting in much larger errors in the derived 
spectral slope of forward-scattering particles for absorption-only models. Similarly, 




16 



fast changes of the single- scattering albedo with wavelength can also cause signif- 
icant errors in the derived spectral slope of cloud and haze particles. 



We also compared our Monte Carlo calculations for three-dimensional, spherical 
planetary atmospheres to calculations for atmospheres composed of plane-parallel 
slabs, with the same optical thickness as the slant light paths through the spherical 
atmosphere, and with scattering pr operties weigh t ed by t he slant op t ical t hickness 
of each atmospheric layer (see e.g. Hubbard et al.l (l2001|) : lKim et al.l (|201 If) for ex- 
amples). For the cases we considered, we found these two calculations to match 
within 10% for optical thicknesses less than 10. This has the implication that the 
contribution of forward- scattered light to transmitted fluxes can be fairly well mod- 
elled using a plane-parallel assumption, for which radiative transfer routines are 
available that are much faster than a three-dimensional Monte Carlo code. 
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